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Abstract 

We discuss problems of functional integral formalisms in a constrained 
fermionic Fock space. A functional integral is set up for the Hubbard model 
using generalized coherent states which lie either in the constrained or in the 
full Fock space. The projection for the latter states is implemented through 
a reduction of the charge fluctuations which induce transitions between the 
constrained and full space. The Lagrangian is expressed in terms of two com- 
plex fields representing spin and charge excitations, and one Grassmann field 
signifying hole excitations. Here, the charge excitations denote transitions be- 
tween states with empty and doubly occupied sites. The projection method is 
inspired by the observation that the local interaction in the model resembles 
a magnetic field in the space of charge fluctuations. Hence the projection is 
understood as an inflnite magnetic field in a spin path integral. 
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I. INTRODUCTION 



Models of interacting fermions on the lattice like the Hubbard model and its derivatives, 
the t — J and Heisenberg models, have been under intense investigation in recent years [|l| . 
These models are thought to describe much of the qualitative physics of the Mott-Hubbard 
metal-insulator transition, including the magnetic correlations and possibly superconductiv- 
ity in its vicinity. Most importantly, there is reason to expect new types of ground states 
such as spin liquid states in the two-dimensional systems. Progress in assessing the valid- 
ity and relevance to many existing ideas proposed in this context has been slow due to the 
lack of reliable theoretical methods. The generic problem here is the dynamics of fermions 
in constrained Fock space, excluding double occupancy of lattice sites. A convenient rep- 
resentation used in this case is in terms of slave particles |Q, assigning fermion and boson 
creation operators to the various possible states at a given lattice site. The occupation num- 
bers of the slave particles at any given site then have to add up to one for all times. Local 
constraints of this type are difficult to handle. Most treatments use a mean field description 
by which the constraint is only satisfied on the average [§,0. Additionally, the slave boson 
methods exhibit a local gauge invariance which is destroyed in the mean field description. 
The conserving slave boson approach presented in [0, however, preserves the local gauge 
invariance. 



One strategy to avoid slave particles builds on a Hubbard Stratonovich decoupling of the 
interaction in terms of density and spin excitations. But a straightforward saddle point eval- 
uation violates certain symmetries, e.g. the spin rotation invariance. Schulz [0 considered a 
time dependent quantization axis for each spin which is being treated as an additional set of 
angular field variables. This procedure preserves the spin rotation invariance of the saddle 
point Lagrangian. 

A different strategy is to try to work in the constrained subspace only. Conventional many 
body formulations build on equation of motion and diagrammatic schemes P-p!0[|, and often 
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involve uncontrolled factorizations which can violate sum rules and conservation laws (except 
for certain limiting C3iSGS, clS the limit of large spin degeneracy in |T0| ) . The use of resolvent 



operator methods may be well adapted to certain problems. Also Quantum Monte-Carlo 
techniques, which work in the constrained subspace only, were devised in |T2|,|13[. 

In this paper we explore the following alternative: the construction of a path integral 



formalism for Hubbard operators ||14| , p!5[| , different from the usual fermionic representation 
||T6| . The path integral formalism has several well-known properties: 

(i) it allows a time-ordered perturbation theory to be developed 

(ii) collective behavior may be incorporated from the outset 

(iii) the classical limit may be easily extracted, e.g. in the case of large spin degeneracy. 

In addition the path integral we develop here does not exhibit the local gauge invariance 
generated in slave boson approaches. 

We organize the paper as follows: 
First, we introduce functional integral formulations in the constrained Fock space. We 
note some problems in the literature concerning the measure. They may be traced back 
to difficulties with the resolution of the unity with normalized coherent states. Therefore 
we review the definition of coherent states in order to extend it to graded algebras. We 
derive coherent states in the full Fock space using two complex fields in order to describe 
the spin and charge degrees of freedom, and one Grassmann field in order to implement 
fermionic statistics. Working with these coherent states, we are able to resolve the unity 
and set up a functional integral with the Trotter formula for which the measure is known 
exactly. In order to test its validity we calculate the so-called atomic limit. To handle the 
problem of a constrained space we develop a well-defined projection method, first for spin 



1/2 functional integrals [p!7| , p!8| in strong magnetic fields, and later adapting it to the electron 
path integral for strong local interaction. We investigate the Hubbard model in the limit of 
strong Coulomb interaction to derive an effective Lagrangian. 
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II. FUNCTIONAL INTEGRAL FORMULATIONS OF THE HUBBARD MODEL 



As a model system we consider a fermionic lattice Hamiltonian with local interactions which 



is most suitable for a discussion of projection methods. The Hubbard Hamiltonian |g;19 
allows "to turn on" the projection by sending the interaction parameter to infinity. Using 
canonical electron creation and annihilation operators it reads: 

<^={T,i} (ij) i 

where is a summation over nearest neighbor sites, and we will assume that tij = t for 
all nearest neighbors on a hypercubic lattice. 

This model has a graded structure in the following sense: 
local excitations of the system either belong to a bosonic or fermionic type. As bosonic 
excitations we classify all excitations involving an even number of electron creation and 
annihilation operators. Spin, particle hole and pair excitations are in this class and will be 
attributed to complex fields in the path integral. Fermionic excitations, with an odd number 
of electron operators, are creation or annihilation of electrons at empty, singly or doubly 
occupied sites, and will be taken care of by Grassmann fields. 

To set up this formalism on the Hamiltonian level first, we introduce local operators 
which represent these excitations more faithfully. Knowing the algebra of these operators 
(the "Hubbard operators" ), one is in the position to derive coherent states for the excitations. 
These Hubbard operators are defined as projections: 

:= (1) 

where \i; fi) are orthonormal states at site i, representing the empty (/i = 0), singly occupied 
(/i =t, i) and doubly ( /i = 2 ) occupied sites (we chose the convention |2) = C|cj|0)). 
Then one can verify the following supercommutator relation: 



s.Axfd^^-i-if'^fxrSpA (2) 
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where S = 2 (&{x^'^ + x"^ ~ 3/2) — 1/2^, and the graded characters x'^'^ of the Hubbard 
operators are in the case of the bosonic operators X^^, X""""", X*^^, X^'') and 

1 in the case of the fermionic operators X^'^, X'^'^, X^'^, X'^^. The algebra of the Hubbard 
operators has been extensively studied in 0. 

Now we may write the canonical creation and annihilation operators in terms of linear 
combinations of Hubbard operators: 

d = X"° + aX^-'' (3) 
However, this transformation is not linear and the inverse of it reads: 



X<rO 


= (1 - n^.) 


XOa 


= (1 - 






X2a 






= (1 - n_a-) 


Xc^-a 


= ctc^. 


XOO 


= (1 - n^) (1 - n_^) 


X22 


= n^n_„ 


X02 


= C^Ci 


X20 


— /^t/^t 



(4) 



The bosonic number operators automatically satisfy the local completeness relation U = 
XOO _j_ j^TT _|_ ^ii _|_ j^22 ^j-^icii follows from the orthogonality and completeness of the local 
states. 

Furthermore, the Hubbard Hamiltonian takes the form: 

H = -tY: E(^r° + ^xtmx^" + ^xr'^') - /i E E + (t^ - 2/i) e xF (s) 

<ij> f i cr i 

Obviously, the kinetic energy is split into four terms, each of them representing a specific 
interaction process between excitations on nearest neighbor sites. The local terms in this 
Hamiltonian may be interpreted as chemical potentials fi and 2/i — f/ for singly and doubly 
occupied sites, respectively. 

Now it seems very easy to project out doubly occupied states with these Hubbard opera- 
tors . Since the last term in @ attributes an exponentially small weight to doubly occupied 



sites for large U it appears to be obvious that the infinite U hmit is accomphshed by discard- 
ing all terms with index 2 in equation We will show that in the Lagrangian formalism, 
i.e. in the path integral formalism, such a projection scheme fails - except for very exotic 
choices of the measure of the path integration. 



Construction of a path integral for the degrees of freedom represented by the Hub- 
bard operators requires a mapping onto classical fields or Grassmann fields depending on 



the statistics. Wiegmann formulated the path integral representation using geometric 
quantization. First, he introduced the projection operator Q onto coherent states. The 
Lagrangian is: 



L = -Tt{HQ) + \ C du Tr(Q duQ drQ) 

Z JO 



(6) 



where Q(t, u) is an arbitrary differentiable extension of Q{t) — only restricted by the bound- 
ary conditions Q{t,0) = const and Q{t,1) = Q{t) and t,u are imaginary time variables. 
The second term in (^) is the so-called Berry phase. This phase is being collected by the 
system while transporting it along a considered path |2^. The first part, containing the 
Hamiltonian of (|^), is the 'Boltzmann weight' of this path. 

In order to investigate the limit of infinitely strong local interaction {U = oo) Wiegmann 
proposed the following form for Q{t): 



Q(r) = (1 - V) E + E + E ^°''4v^ + ^ 



(7) 



where 



I 



cos('i9) 



and A„„i 



I 



sin(?9) cos(t9)e^'^ ^ 



V 



sin^('(9) 

sin(^9) cos(?9)e-*'^ co^^{d) 

and V'^V' cire Grassmann numbers. Only empty and singly occupied states are involved in 
(0), making it plausible that a projection onto a constrained space has been successfully 
implemented. Using (P) and (0), the partition function is: 
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Z = Tre-^" = jv 
X exp 







(8) 



(-20 sin^ ^ + [dr + 20 sin^ + Tt{QH)^ dr 



To point out the problem with the normahzabihty of Z we consider the following trivial 
limit Tt{QH) =0. In this case the partition function should be equal to the dimension of 
the Fock space. Applying the relations of appendix |B|, it is straightforward to integrate out 
the Grassmann degrees of freedom, and we find: 



I?[^(r),y,(r)] 









1 + exp 


/ dr iip sin 






Jo 





lim (2^ + 2) 



(9) 



using the results of appendix where N is the number of time slices in the discretized 
partition function. This partition function cannot be normalized to the correct result, 1 + 2, 
with a common prefactor in the spirit of Feynman's path integral. 

To further illuminate this problem, we investigate the measure. It should be possible to 
resolve the unity operator with coherent states: 



^ d§ d^ I dip I dif)^ ^, ip, ij^) g = + + X^^ = 1, 



(10) 



where /i is a measure to be determined. It can only be a function of the spin variables, 
otherwise the term involving X^^ would be eliminated. From (|^ the prefactor of in 



d^ ("^^ d^ /i(^,<^) = 1 
Jo 



(11) 



and accordingly of X^\X^^: 



r27V rn rlir 

d§ / d^ /i(79, (/?) sin2(^9/2) = 1 d§ d^ /i(i9, y;) cos^{d/2) = 1, 

Jo Jo Jo 



(12) 



Adding the two last equations results in a contradiction to (|TT]). As a consequence, the 
resolution of the unity operator is not possible with Q of equation (0) []. However, the 



-^It seems, the doubly occupied sites have not been projected out consistently. 
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derivation of (H) requires such a resolution. 



Weller ||T5[ was aware of this problem and succeeded in setting up a resolution of the 
unity operator. He obtained the following relation: 



Q 



1 + + i) sin(t9/2)e^'^X^° + cos(^9/2)X^° 



|0) 



1 + -ij^j^ + X^^^ sin(79/2)e-*'^ + X^^^j cos(^9/2) 



^00 ^ ^TT + x^i 



1 



(13) 



where 



(14) 



But Q of |TH] is not a projector, i. e. ^ Q. In other words the coherent state is not 
normalized. Using the relation (0), Weller derived a path integral. The resulting Berry 
phase has an unusual form: 

1 TV N 



^ n=l n=l 



X 



X (sin(?9„/2) sin(?9„„i/2)e*('^"-'^"-i) + cos(^9„/2) cos(i9„_i/2) 
/ dr (d^ - ,0 sin2(^/2)) ^ - - E ^i^n 

"'O n=l 



(15) 



It is impossible to take the continuum limit because the last term is missing a prefactor 
e = (ijX . The Hamiltonian part, Tr{HQ), has the same form as in Wiegmann's paper. 
One can prove that a resolution of the unity operator with normalized coherent states is not 
possible in the case of the U = oo algebra (see appendix 



To investigate the U = oo limit, Schmeltzer |£T[ used a slave boson technique in the full 
Fock space: 



L — Lq + Lx + Lh 



(16) 



with 



^o = E 



Lh = y: - - E E [uM^lh.i^i + H.c.) (17) 

i o" <i,j> 

where are complex fields representing the spin degrees of freedom, A is a Lagrange mul- 
tiplier field, and ip is a Grassmann field for the hole excitation. 
A graded transformation was set up 



= z^^l - A = A (l - ^V) (18) 

which transforms the graded constraint into a trivial complex constraint which can be solved 
by parameterization. Consequently, the Lagrangian is L = Lq + -La + ^-h^ where 

i 

Lx = ^iA, (^zl^Zi^^ - 1^ 

= E [f^'^i'^i - + E {tiMzjaZja-ipj + H.C.) (19) 
i <i,j> 

The resulting Lagrangian is identical to Wiegmann's, and therefore suffers from the same 
problem of normalizability. 

Schulz also worked out a path integral formalism for the full Fock space. He started 
with the canonical formalism and decoupled the model by a Hubbard Stratonovich trans- 
formation. Having introduced a unity at each space and time step, he chose the unity to be 
1 = RR\ where is a spin rotation matrix SU{2)/U{1). In order to make R a dynamical 
field, he additionally integrated over all angle variables at each space and time step using 
the invariant measure of the group. In this way he built in a dynamical, rotating spin ref- 
erence system. After taking a saddle point approximation with second order fluctuations, 
he applied a 1/U expansion and integrated out the Hubbard Stratonovich field variables. 
Again the effective Lagrangian determined in this way is equivalent to Wiegmann's in the 
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limit U = oo. Also Weng et al. introduced a dynamical rotating spin reference system 
in the same manner as Schulz. 



III. COHERENT STATES 

In the following sections we will derive a path integral formalism which uses only one 
Grassmann field variable and works with a well-defined measure. It should give a correct 
result for the partition function and, in addition, the measure of the path integral should 
possess a controlled continuum limit in the time direction. 

First, we reconsider how to set up the definition of coherent states. A well known 
procedure is to require the coherent state to be an eigenstate of an annihilation operator. 



This results in an unnormalized state (see ||23|). A straightforward generalization of this 



scheme to an arbitrary operator algebra is unclear. Therefore we take the definition of 
Perelomov p^j25|| for generalized coherent states. The underlying idea is that the coherent 
states are the orbit of a Lie group represented in a Hilbert space and acting on a reference 
state. This scheme produces normalized coherent states. We will discuss only those coherent 
states which fulfill the two conditions: (I) normalization and (II) resolution of the unity 
operator. This will automatically result in a well-defined Berry phase in the continuum 
limitQ. 

The limit of infinite local interaction is of special interest in this paper. It should be 
advantageous to construct normalized coherent states in the constrained space. However, 
we will show in appendix |C 2| that possible coherent states with one Grassmann field in the 



constrained space will not fulfill (I) and (II) simultaneously. 

For the coherent state to be defined uniquely, up to a phase factor, the Lie group should 
be divided in left cosets of the maximal isotropic subgroup, i.e. all operators, which when 



^Abandoning condition (I) may result in a Berry phase which has no continuum Umit, cf. Weller 
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applied to the reference state generate just a phase factor, are removed. 
In the case of bosons the coherent state is: 



gaat-aat|Q^ (20) 



and in the fermionic case: 



e'^^'-^^'lO) (21) 

where a is a complex variable, ip a Grassmann number, a, bosonic creation and annihila- 
tion operators and C, the fermionic operators, respectively. The unity operator U does 
not appear in the exponentials of equations ( ]20|) and because it is a generator of the 
maxumal isotropic subgroup. 

The graded case is more involved. As mentioned above, bosonic and fermionic Hubbard 
operators are the generators of the graded Lie group. To get the new coherent state \G) we 
write the analogous ansatz: 



\G) = exp 



T) (22) 



where we have just chosen | f) to be the arbitrary reference stateQ. Now, we intend to use 
only one Grassmann field and the smallest number of complex fields possible. The minimal 
number of complex fields is two. This expresses the fact that we deal with a spin and a charge 
degree of freedom. Therefore, we choose xto = c^^, Xot = Xio = Pi^^ Xoi = /^V^ 

and a, (3, Xab to be complex. Then the exponential function is expanded using ip"^ = 0, etc. : 

\G)= U + T.XabXab+lY.Y.XabXcdXabX,, + ...] |T) 
y ab a^b c^d j 

= (a + I T) + (c + li) + eV^ |2) + |0) (23) 



^The discussion of Klauder |17| concerning the correct choice of the reference state is not applicable 
here because both | |) and | |) are states with maximal weight. 
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where a, b, c, d, e, f are complex variables which are functions of the XaftQ- The conditions to 
be fulfilled by the coherent state are (/) normalization and (//) existence of a resolution of 
the unity operator: 

(/) {G\G) = a^a + a^bip-^^ + afet^^t + c+c + crf^^t + ctrfV-^t + g^e^t^ + /V^V = 1 

(//) 1 = Jdij{a,b,c,d,eJ)dijUilj\G){G\ 

= Jdfi (a, 6, c, d, e, /) ((6at + ab^) X^^ + \e\^X^^ + l/pX^o + (ccit + dc^) X^^ 
= + X^T + + = 1 

Here, fi is an unknown measure which still has to be determined: The prefactors of the 
diagonal Hubbard operators are integrals and should equal 1: 

/ dfi (a, b, c, d, e, f) (ba^ + ab^^ =1 J dfi (a, b, c, d, e, f) |ep = 1 

, / \ - ^^^^ 

/ dfi (a, b, c, d, e, /) |/p = 1 / dfi (a, b, c, d, e, /) icd'^ + den = 1 

In addition, the normalization condition (/) has to be fulfilled. The solution of this system 
of equations is not unique since there are six unknown variables and one unknown measure 
function /i. One solution is, e.g. : 

\G) = (l + 1^^"^) sm{^)e'^ I t) 
+ (l + i^^t)cos(^9) li) 

(25) 

+ Vcos(^) |0) 
+ V^sin(^)e^<^ 1 2) 

This state is normalized and one can resolve the unity operator and obtain a trace: 

11 = 1 \G){Gl Tt{A) = J {xG\A\G) (26) 



^In appendix ^ we will show the calculation of these functions explicitly in two special cases. 

^ The integration over the phases of and {ap + be^) which are the prefactors of the offdiagonal 
Hubbard operators are supposed to give 0. 
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where / = 2 Jq ^ J^'^ ^ Jq ^ J^^ ^ / dil)~^ dip . x indicates that Grassmann numbers will 
acquire a minus sign. In this solution we have used a so-called reduced measure (see 
||15| , |26|| ) . There is another resolution and a trace with the invariant measure / = 
2^ Iq^"^ dQ^m.{2Q) Jq^ d(f) Jq^'^ d')!}sm{2')!}) Jq^ d(p J dip^di/j; however the Lagrangian will be in- 
dependent of this choice. 

The expectation values of the Hubbard operators are listed in the following table. 



bosonic 


fermionic 


(G|XOO|G) = 


= ip'^ip cos^ (6) 




(G|XTO|G) 


= ^ sin (^9) e-^^ cos (0) 


(GIX^T \G) 


= (l + ipip^) sin^ (d) 




(G|XOT|G) 


= sin [d) e'^ cos (^) 


{G\X^^ \G) 


= (l + ipip^) cos^ {^) 




(G|XiO|G) 


= ip cos (i?) COS {6) 


{G\X^^\G) -- 


= ijj^ip sin^ {6) 




(GlXOilO 


= ip'^ cos cos {9) 


(G|X02|G) = 


= ifj^ip sin (6) cos (6) e'^ 




(G|XT2|^) 


= ifj sin {^) e-'^ sin (^) e*"^ 


(GlX^oiG) = 


= ip'^ip sin (6) cos (6) e~^'^ 




(GIX^TIG) 


= ip^ sin (^9) e*^ sin (6) e"^"^ 


(GlX^i \G) 


= (l + ipip"^^ sin (-i?) cos 




(G|Xi2|^) 


= cos {'&) sin (^) e*"^ 


(G|XiT \G) 


= (l + ipip^^ sin cos 


e*^ 


(G|X2i|G) 


= ip^ cos (^9) sin (6) e'^^ 



This table allows y?) to be interpreted as angular spin variables, and [6, cp) as angular 
charge variables. The charge variables are pseudospin variables in the space of empty and 
doubly occupied states. They parameterize rotations in this space, as (^9, ip) do in the spin 
space. The Grassmann field induces transitions between spin and pseudospin space, e.g. 
creates a hole in the constrained space and a doubly occupied state in the complementary 
space. 

IV. PATH INTEGRAL REPRESENTATION 

In the previous section we introduced coherent states for the Hubbard operators and fixed 
the measure which guarantees a resolution of the unity. Now, using the Trotter formula and 
this resolution of the unity operator with the coherent state ), we are in the position to 
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set up a path integral for the partition function. The derivation is analogous to the spin 1/2 
case in appendix 

r/3 



exp 









where the Berry phase Lq is: 

Lo = {G{T)\dr\G{r)) = [i^, sin^ + 4 {d^ + sin^ 6, - sin^ ^, 



(27) 



(28) 



The Hamiltonian part may be directly taken from the table. It contains the Coulomb 
potential UX'^'^ and the chemical potential \i (Z^a-^'^'^ + 2X^^) which corresponds to 



(29) 



The form of Li allows to interpret — f//2 as a magnetic field coupling to the pseudospin 
(c.f. appendix 0). The discrete version of the local Lagrangian Lq + Li, will be presented 
in the next section. 

Further it contains the kinetic term which we write in the discrete version in order to 
exhibit the correct time ordering of the field variables: 



Lt — Lt-oi + Lt-20 + Lf^ii + Lf^ 



21 



(30) 



where 

Lffil= Uji^i,n-li^j,n cos 9i^n-l cos 9 j^n 

<i,j> 

X (^sin'(9j_„sin^9j^„_ie"'^'^"'""'^-''"-i^ + cos^9i,„ cos^9j,„_i^ 
represents the exchange of a singly occupied site with its empty nearest neighbor, 

Lt;20= Y UjA,n-li'j,n~lCOs9i.a_lSm9j^n-ie"^''"~^ 
<hj> 

X (sin di^n COS 'i9j,„e~*'^''" — cos di^n sin ^j^n^^"^^'' 



(31) 



(32) 



signifies a transition from a state with two singly occupied nearest neighbor sites to a state 
with neighboring doubly occupied and empty sites, 
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-i4>i,n 



Lt-n= ^ii^|,n^],n sin 6'j,„ cos „e" 

<id> 



X 



is the inverse process, and 

Lt-2i= J2 ^^iV^^>V'J>-lSin^i,„sin^J•„_le-^('^''"-'^^'"-^) 

<i,j> 

X (cos^9,,„_i cos?9j,„ + sin?9i,„_i sin^j^ne''''^''"''"^''"^) 



(33) 



(34) 



denotes the exchange of a singly occupied site with its doubly occupied nearest neighbor. We 
check that, if we set all angle- variable fields to zero, we find a Lagrangian which describes 
spinless fermions up to a particle hole transformation (cf. appendix P 



V. ATOMIC LIMIT 



A. Continuum Version 



The partition function may be calculated easily in the case of no band motion: tij = 0. 
Although this atomic limit is trivial physicswise, we consider it carefully in order to test the 
integration techniques and further to develop projection methods in the next two sections. 
In this atomic limit we integrate out the Grassmann field using appendix |^. 



Z. 



at 



V 



exp 


l-i 







- / dr (Lo + Li) 



P[0(r),0(r)]exp 



- J^^ dr (^z0sin2 e+(^fi-^) (cos(20) - 1] 



(35) 

P[t9(r),^(r)] 



The first term is a path integral for the free spin 1/2 and, therefore, yields 2e^^. The second 
term is also a spin path integral, but with the pseudo angle variables and subject to a 
'magnetic field' ^ — U/2. Therefore we obtain 



Z = 2e^^+fl+e(2'^-^)^ 



(36) 
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B. Discrete Version 



As is well known, the continuum version of the action is only a short form, correct to first 
order in (3/N, but several situations necessitate to consider the second order, e.g. magnetic 
field problems or particle hole transformations (appendix [B 11 ). Therefore we reformulate 
the previous subsection in the discrete language and introduce some notation which will be 
used in the following. 

The vector notation of spin coherent states at time step / is: 



|n.) = 

and for the pseudo spin coherent state: 

|N,) = 



(37) 



sm(^y,je 

To first order in e = 13/N we are allowed to use: 



(38) 



1 — (ni|ni_i) = iip sin^ 'de + C(e) 
1- (Ni|Ni_i) =i0sin2^e + C(e) 



(J<ii\a,\Ni_,) = -cos{29i) + 0{e) 



(39) 

(40) 



With these relations we rewrite the atomic limit in the following form: 

N \ N 



■'at 



(41) 



\i=i 



1=1 



X exp 



N 



J2 [^li^i + i^Hi-i (-1 + (nz|nz_i) - (NzlU - ^le + (/i - |)a,e|N,_i) 



Expanding the exponential function and exploiting the fact that all higher orders vanish 
because of the anticommuting properties of the Grassmann fields ipiip'^ leads toR: 



^In functional integrals over Grassmann fields it is necessary to include the zeroth time slice 
explicitly due to the antiperiodic boundary conditions. However, in order to present the following 
calculations in a more readable form we will always skip this term. 
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N / / U U 

X n + ^i^i + ^i^i-i - (n/|ni_i) + (N,|U - -He + (/i - -)a,e|Ni_i) 



In order to carry out the integrations we first collect the complete Grassmann chains as 
explained in appendix O. 



■'at 



N 



n ^d^idifidOidci^id^id^j ) e^" 



N 



n WV^/(ni|ni-i) 



u=i 



+ n U/^;-i(NHll - -He + (/i - -)a,e|N,_i) 



(42) 



Integration over the Grassmann and spin (first expression) or pseudo spin (second expression) 
variables yields: 



1 \ N 

\i=i ^ / 1=1 

r f ^ 1 \ ^ U 

+ / n —dOid<pi e^'^ n(Nda - irle + 
•' \i=i ^ / z=i ^ 



(/i - ^)(7,e|Nz_i) 



(43) 



Using the resolution of the unity operator we immediately obtain: 



Z„t = Tr(ll)e^^ + e'^^Tr exp 



-|ll/?+(/i-|)a./? 



(44) 



VI. METHODS OF PROJECTION USING SPIN COHERENT STATES 

In this section we introduce several methods of projection with the intention of applying 
one of these to the path integral (PTj). It was shown in section [V] that the local Coulomb 
interaction U is similar to a magnetic field in the pseudo spin space. Therefore, we first deal 
with the system of a 'spin in a magnetic field' and investigate the limit of an infinitely strong 
magnetic field. This will result in a procedure with a well-defined measure which allows to 
take the limit of strong interaction U — *■ cxo. 
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A. Method 1: Direct Integration 



We write the spin path integral in its discrete version 



hm (\ 



N 



Y[ dipnd^n exp 



Vn=l 



N 



n=l 



(45) 



and carry out the angular integrations. To this end, the ip integration is extended to L 
intervals and sin^('(9„) = | (1 — cos(2't9„)) is replaced. We find: 



Z= lim lim — n / 



2nL 2L Jo 



d'dr, 



X exp 



N 



n=l 



(46) 



where periodic boundary conditions were used. Integrating out the N {} variables, a A^-fold 
convolution product of Bessel functions Jo emerges which may be diagonalized by Fourier 
transformation. Furthermore we integrate out the ipn and use the resulting 6 functions. It 
follows: 



/•1/2 dk 
lim c ■ / — 

N-*oo J-1/2 277 



1 



N 



(47) 



where c = (l/vr)^. Taking the limit — > oo, the only surviving contributions are due to 
both poles at k = ±1/2. In order to determine the contributions, we introduce a cutoff ec- 
Integration over k yields: 



2 lim lim 



(48) 



AT^ooEc^o ec7r(A^ + 2) 

The prefactor 2 arises from the two poles. To solve the problem of the continuum limit in 
the spirit of Feynman requires only the quotient of path integrals with different couplings 
to be well-defined. The existence of this quotient is guaranteed by the precise knowledge of 
the normalization factor. In our case the limiting procedure is explicitly known (see (^Sf)). 
Specifically, the problem of a spin in a 'magnetic field' is solved analogously to (|45|)-(^8D 
and yields the same normalization factor as in ( P8| ) — independent of the 'magnetic field'. 
Therefore the existence of quotients such as Z{B)/Z{B = 0) is ensured. 



B. Method 2: Time Dependent Saddle Point Approximation 



The spin path integral is the sum over all paths on a sphere. A path is a set containing 
points for each time and space step. In order to map the surface of the sphere onto planar 
coordinates uniquely the sphere is cut along the equator and each of the hemispheres is 
represented by a separate plane. With this stereographic projection each point will be 
placed on either of the two planes. 

There is a natural way to extend the standard saddle point approximation to this spin 
path integral: at every space and time step the functional will be expanded around a point 
in each of the two planes. In this section we will clarify this procedure and show that it 
results in the exact expression for the partition function. 



For the explicit evaluation of the partition function of a single spin, first, we chose north 
and south pole as reference points for the expansion: 



Z — lim 



N-*oo VtT 



N 



N 



Yl dipnd^n 1 exp 

n=l / 
^ ' n=l n=l 



N 



X) ( (<^n - <^n-l) 2 COs(2i?„) 



n=l 

-i(ipn-<Pn-l)(h-^l 



(49) 



+ e 



+i(<^„-<^„-l)(i-1?2 



where we have expanded the exponentials up to second order in where the field is the 
fluctuation around north and south pole. Next we integrate out the The integration runs 
from to arbitrary p. The result is an infinite convolution of Presnel functions which can be 
diagonalized by a Fourier transformation. Making use of the periodic boundary conditions 
for the partition function reads: 

1 1 



lim c- dk 

N-^oo J 



+ 



N 



lk+\ ^k-\ 

The only contributions to the integral are due to the poles at ±|. 



(50) 



One might worry that c has to be tuned separately for each kind of spin system. Hence 
one might infer that the quotient of two path integrals is not well-defined. Therefore we cal- 
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culate two spins in a magnetic field interacting via a Heisenberg coupling. Using the above 
saddle point approximation and introducing auxiliary fields, to write cos('i?„) as a derivative 
to tpn — ^n-i, we show that the same limiting procedure, using c and cutoff ec, normalizes 
this more complicated problem. 

First, we present the calculation for an Ising coupling of two spins in the presence of a 
magnetic field. The Lagrangian is: 

IT' J 

L^-T. ^ cos(2^?„,,) + - cos(2^?„,i) cos(2^?„,2) + B J2 cos(2^?„,0 (51) 

i=l,2 ^ i=l,2 

where Xn,i — (fin—fn-i- Now we work with the Gaussian method discussed above. Therefore, 
we have to add up the four combinations of north and south pole at each site. Rewriting 
the Ising coupling and magnetic field with auxiliary fields, the partition function reads: 

2N N ( 

n n d^n,idl^n,i \ Yl 12 

I <Tl=±lCT2=±l 



1=1,2 n=l 



X exp 
X exp 



(52) 



Next we have to do the i? integrations which cover the interval from to p. Fourier trans- 
formation and straightforward calculation yields: 



Z — lim ■ 

AT— ►oo 



dki / dk2{f{ki,k2)) 



N 



(53) 



where 



f{ki, k2) = ^Yl [e""^'^^ + e-"'-^'/^ (cos(iJe/2) - 2icTi sin(iJe/2)) 



X 



-2Be 



X- 



(54) 



^1/2 + a2ki^l/2 + (7i(72A;2 
The dependence on p has vanished because we have applied formulas of the following kind: 



/ 



°°dxSiJx) .b'^x. Jt^I 
2p ^ "2 ^ 



-5/2 



P 
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(55) 



Using the same cutoff and normalization constant as in method 1, we obtain: 

Z = le-'^l^ cosh{Bf3) + 2e'^'^^^ (56) 



which is the correct result for the model Extending these considerations to a Heisen- 

berg coupling, we switch on the x — y interaction, expand in the x — y coupling to all orders, 
calculate the value of an arbitrary term, and sum up all terms to leading order in e. Using 
the same limiting procedure and normalization, we find the correct result for the partition 
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Later on, we will need a well-defined projection method to handle the U ^ oo limit. 
Thereby we note that the north pole collects a weight e^^ and the south pole e~^^ . The 
only surviving path in the limit of an infinitely strong magnetic field is the path that stays 
at the north pole over the full time range. However, if we set "i? = 0, we would not find 
the same normalization constant as before. The correct way of doing the projection is the 
following: 



L = i(p sin^(?9) + B cos{2§) L = i(p sin^(t9) 



+ B 



(57) 



north pole 

In this prescription the Berry phase is not eliminated, but has some small fluctuations around 
the pole. In the Hamiltonian part, the prescription yields the same result, as if we had set 
d = without further considerations. We note that this prescription does not change the 
measure. 

In the discrete language the corresponding prescription is: 

\i=i ^ / 1=1 

/TV ^ \ N 

~^ / n —d^id^pi m^il (1 + Ba^e) Pnorth poie|ni_i) 
\l=l ^ / 1=1 



(58) 



where Pnorthpole 



. The path integral for infinite 'magnetic field' may be either 



gained by using (^) in the continuum version (analogously to |V A|) or by using (|58|) in the 
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discrete version (analogously to [VB|) . As compared to section 0, the projection onto the 
constrained space is now implemented into the path integral formalism ab initio. 



C. Method 3: Action Angular Variables 



In this section we want to present a method, interesting in a didactical sense. First, we 
consider the Berry phase of the spin 1/2 path integral. 



lim 



N 



Y[ d^nd'^n exp 



\n=l 



N 



n=l 



(59) 



There is a geometrical interpretation to this expression: The Berry phase is the area enclosed 
by a single path. This reminds us of classical action angle variables in classical mechanics. 
Therefore, we want to implement similar variables. 



To do so, we start by separating the partition function into a sum over all right and left 
rotating paths: 



1 /l^^ 
Z = - lim — 

2 N^oo Vvr^ 

+ - lim — 

2 N-*OQ VtT^ 



N 



1 



N 



W di^nd(Pn exp 



Vn=l 



TV 



Yl d'&ndipn exp 



Vn=l 



'i ^ 

^5 XKv^" ~ V5«-i) cos(2^9„) 

.-^ n=l 



(60) 



The classical action variable is the area enclosed by a single path. Therefore the action of the 
single spin is just this classical action variable, and it seems to be advantageous to transform 
the field variables into action variables. Since the field 'd consists of t^ri-variables (and 
there exists only one classical action variable), we introduce new variables being defined as 
the areas which are swept out by the geodesic line from the north pole to the actual position 
on the considered path (see figure). In this way we obtain A^ variables: 



Bri = 7;J2(Vi- Vi-i) cos(2i9/) 
^ 1=1 



(61) 



Then is the Berry phase. The Jacobian is: 
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IJI 



dB 


dB 




dB 


dB 


d-d 


dip 




d^ 


dip 


dip 


dip 







1 




dip 





(62) 



where 



d'drr, 



{(fim-i - Vm) sin(2t9„) m <n 
m > n 



(63) 



is a triangular matrix and we find explicitly: 



N 



N 



1=1 



1=2 



The partition function reads: 



(64) 



The action has become trivial but the integration over the Jacobian is a A^-fold convolution. 
Fourier transformation decomposes the convolution into a product of independent integrals. 



2 



(65) 



and 



h* = j dk,,N /^(A;i,;v,+1) 



(66) 



where 



f{ki,k2 



xy 



2 I 2 

X +y 



-1/2 



exp [ikix + k2y\ 



rdrdif 



^ikr cos tp 



^ k ^kl + kl ^^^^ 
The normalization of the resulting singularities is similar to the procedure of the section 
above (see E3-HS|). Furthermore, the analogous scheme applied to a spin in a magnetic field 
yields the correct result. 
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VII. ALTERNATIVE PROJECTION METHOD USING FERMIONIC 

COHERENT STATES 



In the first part of tliis section we present an alternative projection metliod using tlie 
standard patli integral formalism with two Grassmann fields. The basic idea is to perform 
projections at each time step explicitly through a suitable choice of the measure. Then, we 
introduce coherent states in the constrained space and show that their measure is identical 
to the 'projecting' measure of the first part. 

First, we repeat the calculation for a local action, e.g. a chemical potential term: 



■'at 



N 



\n=l 
/ N 



N 



X • 



n dip\nd'^]ndipl^dilJin exp 

n=l / 
N \ 
n #Tn^^Tn#In#inj 

N \r/^ \/^ 

n ^T,nV'T,„ n ^lM,n + n ^I,n^i,n-l(l + /^e) 



\n=l 
f N 



\n=l 



\n=l 



\n=l 



N 



N 



[ (68) 
The last 



\n=l / \?i=l 

In the last step of (^) we expanded the exponential as in appendix | 
product of (1^) represents the double occupation (see appendix |B]). Now we introduce a 
prefactor in the measure which destroys double occupation and, simultaneously, allows the 
other three states to survive with the same weight: 

N \ N 



■'at 



n=l 



X exp 



N 



o" n=l 



(69) 



The quartic term annihilates all the terms of the exponential function, except the zeroth 
order term. In order to explore further this projecting measure we re-exponentiate as 



exp [iXn (l + 'ip\,n^hn^in^Ln " ^\,n^hn " ^i,n^i,n)) 



(70) 
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to find a contribution to the action: 



/3 



In the continuum limit e — > 0, diverges. This resembles an infinitely strong local inter- 
action at each time step. After translation to operator language and performing a particle 



hole transformation (see appendix P 1|) one recognizes a local interaction term with infinite 
strength: 

^(U - B^B\ - B^B\ + B^B\B^B\) ^ -^C\C^C\C^ (72) 

This ad hoc scheme of introducing a measure, which projects into a (fermionic) subspace, 
may be embedded in the more general language of section |ITT|: The commutation relations 
of the operators X"'^ (where a,6 = 0, t, |) clearly constitute an algebra in the constrained 
space. So we may set up coherent states as: 



|Go) = exp^T^^°+^i^^°-'^-- |o) 





IT) 






li) 


(73) 



where V'o- are Grassmann numbers. Now the projector onto the normalized coherent states 
reads: 

+Xii(V^iV^I + V^}^tV^iVi) (74) 
For the resolution of unity we have to find a measure ^{ipp ipi^ , if)'^^, if) y) which satisfies: 

di)^di)\d'4)idip\ /i(^}, V^i, t/;}, V^JP = 1 (75) 
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Obviously, any measure consisting of a single product of the Grassmann numbers will be 
inadequate to solve (fTSl). It turns out that ii has to be of an additive form 

V-T, ^I, i^l) = - (i^l^Wl + ^T^T + ^I^i) (76) 

in order to ensure the resolution of the unity operator. 

The Berry phase has to be calculated in the discrete version. Terms in fourth order of the 
Grassmann fields arise, but are eliminated by multiplication with /i of (^). The surviving 
contribution has the canonical form J2c7'4'i9r'4'a, as in (|69D. Having fixed the measure and 
the Berry phase, the functional integral formalism is established. One might worry that the 
derivation of the Hamiltonian terms of the action always involves higher orders in the field 
variables; e.g. the kinetic term as derived with ( |73D is 




Fortunately, the measure /i eliminates all fourth and higher order terms, similarly to the 
calculation of the Berry phase. 

Before we try to derive an effective Lagrangian for large U, we should summarize the 
projection techniques of the last two sections: It turned out that it is possible to construct 
normalized coherent states which lie in the constrained space and resolve the unity. But 
their measure is of an additive form, i.e., a polynomial of the Grassmann variables for up 
and down spin. The minimal number of complex Grassmann fields for such a resolution 
in the constrained space is two, so that the proposed coherent state does not render the 
path integrals more tractable as compared to the canonical formalism with two Grassmann 
fields in the full space: The local interaction is suppressed but at the expense of a most 
inconvenient measure. Therefore, coherent states should be set up in the full space in which 
we may reduce the number of Grassmann fields to one. Additional complex fields have to 
be introduced for spin and pseudospin. Since the local interaction denotes a 'magnetic field' 
in the pseudospin space, we investigated an elaborate saddle point technique for spins in 
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magnetic fields. It serves to constrain the excitations to the considered subspace by sending 
the magnetic field to infinity in a controlled way: The pseudospin fields may be fixed to their 
'subspace direction' everywhere except in the Berry phase where small fiuctuations about 
this direction have to be kept, even for infinite U. 



VIII. CUMULANT EXPANSION FOR LARGE LOCAL INTERACTION 



Now we will show how to gain the effective action to lowest order in t/U using the 



coherent states (pSj) and the projection method of section The effective action will be 
similar, but not identical, to the action derived by Schulz through a Hubbard Stratonovich 
decoupling scheme. 



The Lagrangian for the Hubbard model was introduced in section [V 



L = fi + iipi sin^ i^i + 



+ ipj ['^r + ^ + i4>i sin^ 6i - iipi sin^ + (^/i - ^ j cos(26'i) j ^/'j j + U 
Applying the projection prescription the Lagrangian reads for U = oo: 



U 



U 



(78) 



L = V ( -/i + i0i sin^ + 7/^1 ( 9^ + /X + z0i sin^ 6^^ , , - iifisin^ ^pA + Lt (79) 

^-^ \ \ north pole J J 

where Li restricts hopping to the exchange of a hole with either a |-spin or a |-spin on the 
nearest neighbor site: 



/ ^ Ujipitpj (^sin-i^j sin-i^j-e ^'-'^^ '^•'^ + cos-t^j cos-i^j^ 

<i,i> 



(80) 



The path integral with L from (^) may be cast into a more conventional form without 
a prescription. Through a redefinition of the measure 



V 



^(r),0(r),^(r),^(r),^t(^)^^(^) 



V 



exp 



north pole ; 



(81) 
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the projected Lagrangian acquires the form: 

Lu=oo = Y1 {-f^ + i^i sin^ + tpj [dr + - i^Pi sin^ ipi 

i 

+ ^ tijVi^J [(sin?9iSin?9j-e-*('^'~'^^^ + cos^9iCos^9j)] (82) 

This form of the path integral would be more useful if we were able to integrate over the 
purely kinematical fields Q and first. 



In the next stage, we will explicitly calculate local propagators in order to gain an 
effective Lagrangian which is generated by a hopping expansion. The dynamical processes 
of this effective Lagrangian lie in the constrained space only. We will find the propagator 
by integrating first over the Grassmann fields and then the (pseudo) spin fields. As an 
introductory exercise we calculate Z for the atomic limit with Lq = L — Lt (analogously to 
section [VB|) : 



V 



t9(r),^(r),^^(r),0(r),^(r),^t(^) 

N 





/•/3_ 


exp 


- / Lo dr 




Jo 



/ [\{—ddidipi] X{{ni\ni_^)e^^ + 
\i=i ^ ) 1=1 

+ / n —dOld(t)l e'^^ n(Nd (a + /WO^.e) Pnorth pole|Nz_i) 

\i=i ^ / 1=1 



where Pnorth pole 
chains survive. 



2e^^ + 1 



(83) 



and we have used the fact that only the two complete Grassmann 



Similarly, we may evaluate the local propagators necessary for a cumulant resummation. 
In order to check the consistency of the calculation we compare the propagator obtained 
with L to the projected propagator assigned to L. 



First we expand the functional in terms of the hopping parameter t. Odd orders in 
the hopping expansion vanish because they contain an odd number of Grassmann numbers 
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at each of the sites involved in the hopping process. In second order we average with the 
original action Sat — lo {L — Lt). For example the propagator Z-{T Xj;^Q{rm-i)X^=o{Tn)) 
is: 



V 



\9{T),^{r),9{r),cj>{T),^|;^{T),^|;{T) 

X^m-l^i COs(^^_l) C0S(^„) C0S(1?^) COs(l?„_i) CXp [Sat] 



(84) 



where n < m is assumed. With the identity 

cos(6'n) cos(6'm-i) = (N^-i| 

we get 



^0^ 



|Nn) 



(85) 



\l=m 



a=m+l 



n-1/ 



X 



™-l / /TJ U \ \ 

^ n_ (V'/*-i(N,|Il - (-1 - (/, - -)<T,j £|N,_,>j (N 



|N„_i; 



^0^ 



n— 1 m—1 

UMi n i^i^i-i 

l=m l=n+l 



xTr 



^0^ 



exp 



-^1 + (/^ - -^)(^z ) e(m - n - 1) 



gM/3g-/i(m-n-l)e 



(86) 



where the following convention is introduced: if there is a sum (product), whose upper index 
is smaller than its lower index, the sum (product) denotes a sum (product) from the lower 
index to /3 and from to the upper index. With this convention we can also prove the result 
to be true in the case of n > m. 
By exploiting the identity 
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sin(^y„_ij sin( 



J{<f>m-l-(j)n) 



(N„| 



(io\ 



|N„_i) 



(87) 



we are able to calculate remaining term from the second order cumulant expansion. It which 
will generate virtual double occupations in the effective action below. 



V [&{t), v?(r), e{r), 0(r), ^t(r), ^(r) V-^-iV^^ 
X sin(6'^_i) sin(6'„)e 

(io\ 



"^"^ cos(7?„) cos(?9„_i)e 



n-1 



m— 1 



i=m i=n+l 



xTr 



exp 



-^'^ + - l^)<^z] eijn - n - I) 



(88) 



We have to pay attention to the correct sequence of the limits, i. e. to take first U oo and 
then e ^ (for ji <^U). In the case of large [/, it is an excellent approximation to substitute 
for the last result a Kronecker delta with weight given by the integral of the exponential 
over all positive times r = {m — n — l)t 

1 



U 



3^ ^m—l,n 



(89) 



Now we calculate the analogous averages using the projected action, Sat = So dr (L — Lt). 
The calculation corresponds to the preceding: 



V 



i9(r), v9(r), V^^(r), V^(r)] V^^-i^n cos(^9™) cos(^9„_i)e" 



Sat 



V 



■^(r),^(r),^t(^)^^(^)l ^„_^^teA'/3(n^ 



|n„_i) 



n-l 



n—1 



m—1 



Ili^i^l) n ((n^n^-i)) n ^/V^-i(N/|(]l-/ie)Pnorthpoie|N,_i) 

l=m+l l=n+l 



l=m 



gM/3g-M(m-"-i)<; 



(90) 
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which is obviously equivalent to (|8^ ) -(|86| ) . Therefore we are led to introduce the following 
translation rules regarding the transition from the unconstrained to the constrained (infinite 
U) problem. 

{lpllpm-l COS(0„) COs(^„_i)...)5^, ^ (91) 

{iplipm-i sin(^„) sin(^„_i)...)5^, JJ^n,m-i{--)sat (92) 



The ellipses stand for the various spin variables. (^T]) and (|9^) exemplify how to gain all the 
required propagators. The hopping expansion has to be done to fourth order. After using 
these translation rules we re-exponentiated to obtain the following effective Lagrangian: 

i 

-t J2 An-l'^j,na{n^,nj) 

<ij> 
^2 

+ 77 J2 (V'i.n-l^l.n + '^j,n-l'^j,n) (1 " ^i^j) 
<id> 

+ 77 J2 V^i,n-i^i,n«(ni,nj)a(nj,nfc) (93) 

<i,j,k> 

Here Ln = L — Lt and 



a(ni, lij) = y - (1 + mnj) exp [iA{ni, Uj, z)/2] (94) 

A{ni,nj,z) is the area of the spherical triangle defined by the three unit vectors, and z 
points to the north pole of the unit sphere. < i,j, k > denotes the following convention: j 
has to be a nearest neighbor of i, and k a nearest neighbor of j, excluding k = i. To compare 
this effective Lagrangian with Schulz |]^, we have to perform a particle-hole transformation. 
On the operator level the well known anticommutator relation 

C^C = 1 - CC^ (95) 

holds. In the path integral formalism we work with the corresponding 'anticommutator 
relation': 
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N 



Y[ dipidipj ipiipn-i exp 



\n=l 



N 



1=1 



N 

1=1 



(96) 



(p6|) is proved by using relations of appendix |B 1| . Obviously, the information about the 
anticommutator is contained in the time index (see also lETI). The transformation results 



m: 



C — Co + Ct + Cj + Cpair 



(97) 



where 



^ i 

+?/'| Vi,n-i (/ie - 1 + (ni,„|ni,„_i) - (Ni,„|Ni,„_i) 



(98) 



north pole , 

contains the Berry phase which differs from |]^ by pseudo spin fluctuations about the north 
pole, 



<i,j> 



is the hopping term. 



= Tj H - ^lnA,n - i^in'^j,n) (1 " Hill 
<id> 



(99) 



(100) 



denotes the Heisenberg term, and 

Cpair = Yj V^l,nV'fc,n«(ni,nj)a(nj,nfc) (101) 

<i,j,k> 

is the pair hopping term which represents hopping to the next nearest neighbor site and 
thereby transports a pair of spins. We refer the reader to the literature for a discussion of 
these terms, e.g. 0,0]. 

Schulz derived an additional term: 

1 



4U 



5:(1-V^>,)n? 



(102) 



Such a term is not generated in our formalism because no operation will produce a quadratic 
time derivative. 
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IX. CONCLUSIONS 



In this paper we addressed the question of how to set up a path integral formahsm in a 
constrained Fock space. We only tackled methods which do not involve additional 'graded 
constraints' as in the slave boson method. The projection was implemented through either, 
(i) a suitable modification of the measure of the path integral, or, (ii) through adequately 
generahzed coherent states. 

The first approach allows the use of Grassmann fields which correspond to the canonical 
creation and annihilation operators of electrons. A 'projecting measure' was introduced 
which explicitly projects at every time step. 

Alternatively, one might want to consider generalized coherent states which lie in the 
constrained space, exclusively. Such normalized states were constructed explicitly but they 
resolve the unity operator only if the number of complex Grassmann fields is equal to the 
number of spin states. But even then, the path integration is awkward because the Grass- 
mann valued measure is of an additive form which turns out to be exactly the 'projecting' 
measure of the first approach. This implies that normahzed coherent states should be set up 
in the full Fock space first. We presented such states, involving only one (complex) Grass- 
mann field and the smallest number of complex fields possible, i. e. one for the spin and one 
for the charge degree of freedom. The 'charge field', which here signifies transitions between 
states with empty and doubly occupied sites, was shown to be a pseudo spin coupled to a 
'magnetic field', linear in the local interaction parameter. Accordingly, this approach had to 
establish a procedure how to send the pseudo magnetic field to infinity in order to implement 
the constraint of no double occupancy. This procedure fixes the pseudo spin field at the 
north pole in the Hamiltonian terms of the action but keeps fiuctuations around the north 
pole in the Berry phase. It seems to be a generic feature of such constrained systems that 
forces us to keep the Berry phase fluctuations in the unconstrained space, in order to work 
with 'normalizable' path integrals. Finally, this method was elaborated and tested on the 
Hubbard model to derive an effective Lagrangian for strong local interaction. 
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APPENDIX A: SPIN PATH INTEGRAL 



A spin coherent state is given by [^]: 

/ 

|n) = 



\ 



(Al) 



sin (-(9) e" 
cos (t?) 

Furthermore we can find a resolution of unity and a trace with this overcomplete set of 
coherent states: 



H = J dfi (n) |n) (n| 



1 r'^ p'27T 
— j d-d j d(p 

TT^ Jo Jo 



sin2(^9) sin(i9) cos{^)e''^ 

sin(79) cos('(9)e-''^ cos'^i^) 



(io\ 



^0 1 J 



Tr (A) = ^ r d§ d^ (n|A|n) 
~ JO Jo 



(A2) 



Now a path integral is derived using the Trotter formula and introducing an imaginary time 
index: 

Z = Tr e-'^^ 

= lim Tr (e-'"^le-'"^2e-'" ^N-^e-'^^N-2e-'^^N-le-'" 

. N N 

lim \{d^i{jii) n(nz|l -ei/|nz_i) 
1=1 1=1 

Ylm^j\{d^l{nl) exp - ^ (e(n,|9.|ni) + e(n,|i7|ni) + O (e^)) (A3) 



1=1 



1=1 
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In the last step the continuum hmit was taken. Thereby the discrete time index n; became 
a continuum variable n(r) with t = le and the measure is dn{ni) = ^i^i. The continuum 
path integral is written as: 

r/3 



V (r) , (r)] exp 



(iipsm^ {'d) + (n|i^|n)) dr 



(A4) 



We finish this appendix with a note of warning: Usually a path integral with complex fields 
is evaluated through a saddle point approximation. In the spin case it will break down. 
Consider the spin in a magnetic field proportional to B. The Lagrangian is: 

L = i0 sin^ {§) + B cos(2^9) (A5) 

The general solution of the classical equations of motion restricts the spin dynamics to north 
and south pole if the periodic boundary conditions are taken care of. The classical action 
is Sci = ±-B/?, for north and south pole, respectively. Although the saddle point solution 
yields the correct result, fiuctuations around the classical path are divergent, already in 
lowest order. This originates from an expansion of the action about its maximum at the 
south pole. A path on the sphere, on which the spins live, can be described only by two 



complex numbers at each space and time step. In section methods were devised in order 



to generalize the conventional saddle point approximation for spin path integrals. 



APPENDIX B: FERMIONIC PATH INTEGRAL 

In this appendix we want to derive a path integral for a system with spinless fermions 
using normalized coherent states. We take the following definition for coherent states: 

IG) = e^^'-^'^'lO) (Bl) 

Because of the anticommuting properties of the Grassmann numbers the exponential func- 
tion is trivially determined to be: 

|G)= (l + ^^^t^|0)+^|l) (B2) 
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This state is normalized 

{G\G) = ((0| (l + i^^t j + (i| ^tj (^(^1 + l^^t^ |o) + ^|i) j = 1 (B3) 

and we can easily find a resolution of the unity operator 



= |0)(0| + |1)(1| = I1 (B4) 



and the trace is 
Tr A = 



dijU^ ((oi (i + ^^^t^ - (iiv^tj A ((i + ^v^v^t^ |o) + m) (B5) 

In the same manner as in appendix ^ we derive a path integral using the Trotter formula]]: 



. / N \ r ^ 

J ( n d^idlpn j exp ~ ^"-l) + {Gn\H\Gn-l) 



\n=l 



(B6) 



n=l 

If we only consider a chemical potential fiG'^G, we have to replace (Gnl/ZlGn-i) ^ 
—fiip^ipn-i- The evaluation of the path integral is equivalent to the calculation of a de- 
terminant. Alternatively, we can expand the exponential function. Then we have to collect 
a large number of terms which seem to be difficult to handle. But, only the complete chains 
of Grassmann numbers will survive the Grassmann integrations. Thereby 'complete chain' 
denotes a product of Grassmann numbers over all time steps n^=i i^ni^n- This implies that 
the partition function may be written in the following way: 

z= di^idij^ f n ^ni^i - n ^^n-i (i + 

\n=\ ) \n=\ n=l / 

N 

= 1+ n + = l + e'^^ (B7) 

n=l 

The assumption of the continuum limit to exist requires /u(r) to be a smooth function. 
With this assumption we are able to calculate the partition function of a time dependent 
(chemical) potential: 



^Heed the footnote of section VB 
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N 

n 

n=l 



1 + fin^) = 1 + e 



[B8) 



In order to set up more involved path integrals we would like to investigate the question 
of the minimal number of complex fields. Specifically, one might have the idea to introduce 
an additional complex field in (Pl|). 

= gaVCt-atc^t|Q^ =(^l + laaV^t^ |o) + a^\l) (B9) 
The state is also normalized and we can resolve the unity: 



t= dipUip / dada^^i{a) \G){G\ 



(BIO) 



Here, the measure has to obey the relation: 



dada^ fi (a) a^a = 1 



The following Berry phase arises in the path integral: 



(Bll) 



(B12) 



With our consideration of 'complete chains of Grassmann numbers' we can write the partition 
function: 

^ = / n daldandlpldlpn H «n«nV^n^l " H alttn-lV^lV'n-l (1 + /ic) 
\n=l / \n=l n=l / 

We note that each of the two chains possesses a complete product of the complex fields 
Therefore we can integrate out the a„ and we get the known result 

N \ / N N 



Z = const J i f[ di)\di)n | ( H ^"^1 ~ IT i'li^n-i (1 + Aie) 

\n=l / \n=l n=l 



(B13) 



which can be re-exponentiated: 



V 



^t(r),^(r) 



exp 



— dr {ip^drip — /^V' V) 



(B14) 



In this way we have shown the additional field to be superfluous. 
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1. Particle Hole Transformation 

The particle hole transformation is a special case of an unitary transformation for elec- 
tronic states. On the operator level it reads: 

C ^ B (B15) 

Under this transformation the trivial Hamilton operator H = —fiC'^C transforms to if — *• 
H' = —jiBB^ = —fiH + fiB'^B. So, why not introduce the analogous unitary transformation 
in the path integral formalism 

^n^^i, ^i^^n ? (B16) 

In the path integral fiip^tp transforms into —fiip'^ip. However, we expect the particle hole 
transformed result to be yU — fiip'^ip. Obviously, we have to do further considerations to find 
a Lagrangian analogous to a particle hole transformed Hamiltonian. 



The path integral consists of two complete Grassmann chains. One of them is the particle 
chain HV^/V^i-i and the other is the hole chain The particle hole transformation 

should exchange the two chains. Therefore the particle hole transformation is not unitary, 
but 'almost unitary': 

^ (1 - /xe) , V'n-1 (B17) 

and the Jacobian is|: 

N 

\J\ = \{{^ + ^^^) = e'' (B18) 

n=l 

The Berry phase transforms to (first order in e): 



The Jacobian for the Grassmann numbers is defined inversely in comparison to the complex 
numbers. 
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N 



lim J2 ^ii'^n - ^n-l) 
n=l 

N 

^ lim Y: V^n(l - /ie)(^i+i - ^i) = lim i^lO^n - ^n-i) + Oie' 

N—i-oo \ N^oo 
n=l 

and the partition function accordingly: 



N 



n=l 



(B19) 



N 



W dipldipn exp 



\n=l 



N 



n=l 



. / N \ ^ N 

^ 1^1 y n di^idi^n exp - 5: e (^t (^^^ _ ^^^_^) _ ^^^^t + 

\n=l / L n=l 

If we take the continuum limit, i.e. keep the first order in e only, we find: 



(B20) 



TV 



n=l 



(B21) 



So, transformation ( pi7|) produced the expected e^'^-prefactor. 



APPENDIX C: RESOLUTION OF UNITY FOR INFINITE INTERACTION? 

In this appendix we will investigate the possibility to resolve the unity operator with nor- 
malized coherent states which exist in the constrained Fock space. We remind the reader that 
the resolution of the unity operators with such coherent states would result in a well-defined 
Berry phase in the continuum limit. For this purpose we will explicitly construct coherent 
states with either one or two complex Grassmann fields which are clearly normalized. 

However, it will be proved by contradiction that neither of them will resolve the unity 
operator (further discussions in p6[) - assuming a sufficiently 'simple' measure which is not 
a sum of several Grassmann valued terms. 

1. Coherent state with two complex Grassmann numbers 

We write the normalized coherent state in exponential form as in Perelomov (cf. 
(P^)). The bosonic operators will acquire a complex prefactor, and the fermionic a Grass- 
mann prefactor. In the following formula a is a complex field, but (3 is proportional to a 
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Grassmann number i/j-^ and 7 proportional to a Grassmann number t/'j^. The reference state 
is chosen to be be |0). 

1^) = exp [aX^^ + ^ ^^io _ ^*^iT + ^*xOt + |0) 
= |0)+/3| T)+7l i) 

+ ^{«7l T)-«*/3| + 7*7)10)} 

+ ^ {(-««*/? - /?7*7)l T) + («*«7 - If3*f3)\ i) + (7*«*/? - /5*«7)|0)} 

+ ... (CI) 

The construction scheme is found by expanding to eighth order. The coherent state then is: 
r sin(Hl ^ ^ l-cos{\a\) ^ ^„,.„, sin(|a|) - 1^1 ] | 

+ f sin[H) _ ^ 1 -cos(|a|) ^ ^^^^^ sm{\a\) - \a\ \ 
1 lal a laP I 



sin ( a 


) 




a\ 




a 





sin( 


«l)- 






a 


a* 



- (/^*/^ + 7*7) l^i^ - 2/3W7 '~'"^^'"'^~'"''^n 10) (C2) 

which is proved by induction. This state is normahzed. Next the resolution of the unity 
operator is tackled. Therefore we have to take the projector onto the coherent state and 
integrate with some unknown measure. This measure is either a constant, a single product 
of Grassmann numbers, or a sum of such terms. The single product is ruled out, since it 
destroys pref actors of the diagonal Hubbard operators in the projection operator. The sum 



of products yields the measure of section [VI I| as the only possible choice. Finally, the case 



that the measure is a constant is investigated below. 

The arg(a) phase integration over a is supposed to yield in order to destroy the 
terms with off-diagonal Hubbard operators. The remaining terms, with diagonal Hubbard 
operators, contain integrals as prefactors which should equal 1: 



Xtt : / 2^/?*77*^^^^^^ = 1 (C3) 



sin 


(l«l) 


a 




sin( a ) 




a 




a 


3 
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sin( a ) 


a 




sin( a ) 








a 


3 



cos( a ) - 


-1 + 


a 


2/2 




a 


4 







sm(|a|))^ 




a 


4 



\2n 



+ 2- 



(C4) 



= 1 (C5) 



The square brackets can be written as: 

. , 4 sin(|a|)(|Q;| — sin(|a|)) 



\a\ 



(C6) 



Now, adding the first two equations (|C3| ) , (|C4|) results in a contradiction to equation (pSf). 

To complete the proof one can choose another reference state, for example | T) or | |). 
Carrying through the same steps as before one obtains a similar contradiction. 



2. Coherent state with one complex Grassmann number 

Analogously, we try to find a resolution of unity with one Grassmann number. For this 
purpose we take the calculation from the previous section and substitute = '?/'2- Then the 
coherent state is calculated to be: 

[^ sin(|Q|) . l-cos(|a|) ] 

I l"l "* J 

sin(^_^ l-cos(H) | 

[ |a| a ] 

^sinflal) — lal sinflal) — lal , ^ * , 1 — cosdal) 1 , , 

1 _ 7^-/3 ^' ^, ' ' + ^7 ^ \ ' ' - + 7*7 

a\oL\ \oL\a \ay J 

Again, this state is normalized. To resolve the unity we have to integrate with an unknown 
measure. Using the same procedure as before we get: 

J \aY |a| 



\2 

I 



X"": /2(/3/3- + 77*) ^'"^'"'l''^^ =1 



(CIO) 



Adding the first two lines results in a contradiction to the last. 
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APPENDIX D: SLAVE PARTICLE CONSTRAINTS 



To contrast the projection methods discussed in this paper with slave boson techniques, 
we sketch the treatment of slave boson constraints within a path integral formulation [28,4[. 
The trace in the partition function is constrained to some subspace characterized by g = 0, 



Z = Tt 



e-'^"\q = 0){q = 0\ 



(Dl) 



where q labels the eigenvalues of operator Q. We may write 

dX' 



Z = Tt 



Tr 

27r 



(D2) 



The delta function makes sense if we evaluate the trace in the basis of eigenstates of Q. We 
now assume that Q commutes with H to getj^ 



27r 



-l3H+i\'Q 



(D3) 



The integral runs from — vr to vr because the eigenvalues of the number operator Q are 
supposed to be integers. Now, the coherent state functional integral may be set up in the 
usual way (A = A'//?), 



-Tr 



-f3H{X) 



hm / ir 



1 : e-^^^(") : U : e^^^^^^^ : 1 : e'^^^^^^ : ...U 



(D4) 



where H{X) = H — iXQ. We may expand the exponentials in this Trotter formula only 
because (3X is limited to the interval [— 7r,7r]. It remains to express the unities by a closure 
relation for coherent states. 

We should draw the readers attention to the following two facts concerning the constraint 
term: 



^e.g. a path integral for H = J2<i,j>J2a'^ij^icr^j'^ with constraint X'Q = A'(X]i "-it^^ii — 0) is 
invalid because Q does not commute with H. 
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(i) The constraint is not implemented at each time step separately - nor at the initial 
time step. It was implemented at the operator level and is now 'smeared' over all time steps. 
This contrasts with the projection methods introduced in the main body of the paper. 

(ii) A is a parameter (one for each site), and not a (time-dependent) field. It is not 
correct to make A time dependent (A A„). However, since Q is a number operator, say 
iXQ = iX (j2aO'lt(^a ~ l) ^ "i^J^a ^a^^ctn-i ~ '^^i may perform a gauge transformation 
of the fields $0,^ e^^^^a^ to generate additional terms in the Berry phase: 

= - + (1 - e'(^"-^'^"^) (D5) 

The additional term may be written as i$*$Q,9T-x, if drX stays finite (i.e. if Xn — Xn-i = 
0{At)). This gauge field has to obey the boundary condition x(0) — xiP) = Svrm to keep 
exp(iS') invariant, and A — A + ix(r) =: A(r) (see the discussion in [jl]]). Now we sum all 
paths of A(r) in order to obtain a functional integration of the Lagrange multiplier field 
A(r). Interpretation of this procedure is as follows: The static component of the field A(r) 
enforces the constraint whereas the Fourier components may be gauged away. 
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